Pion and kaon parton distribution functions in a meson cloud model 
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In a model of hadrons in which the hadron wave function at low Q 2 is represented by a hadron- 
like Fock state expansion, we calculate the ir and K parton distribution functions by comparing to 
experimental data. We show that the model gives an accurate description of available data and, in 
addition, gives a clear prediction of the initial flavors to be considered at the low Q 2 scale where 
perturbative evolution starts 
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The parton content of mesons is not well known due 
to the scarce experimental information as compared to 
the rich and accurate data existing for the proton. The 
meson structure is measured in dilepton Drell-Yan (D- 
Y) production in meson-nucleon interactions. So far, a 
few experiments |2j have measured the valence quark 
distribution in charged pions, while for kaons, the only 
experimental information today available is on the Uk 
distribution in the K~, which has been obtained by the 
NA3 Collaboration @ in the 80's. 

In D-Y dilepton production in the reaction H A +H B — > 
l + l~ , the biggest contribution comes from the subpro- 
cesses q + q — > l + l~, where q, q are quarks in the initial 
hadrons H A , H B . Then the partonic structure of one of 
the initial hadrons can be extracted if the structure of 
the other one is known. Thus, the structure of charged 
pions can be measured in ir^N interactions. Further- 
more, some authors have claimed that adequate linear 
combinations of ir^p and 7r ± n D-Y cross sections can al- 
low the measurement of the valence quark distributions 
in charged pions independently of their sea quark and 
gluon distributions [Q. In a similar way, the u valence 
quark distribution in the K~ can be measured. How- 
ever, the valence s quark distribution cannot be easily 
determined in kaons since it has to anihilate with a s sea 
quark in the nucleon to produce the pair, giving a 
negligible contribution to the D-Y cross section in the 
experimentally accessible kinematic region. 

Then, to know the kaon structure, which is desirable 
since, for instance, kaons are frequently used as beam 
particles in experiments, one has to rely on theoretical 
models to complement the experimental information. On 
this respect, several models have been proposed. Among 



'Electronic address 
t Electronic address 
* Electronic address 



2avila@uniandcs.edu. co 




saiialn uojmiiaiid 


^s.edu. 


CO 


liiagninCicbpf.bi 





them we can mention the model of Ref. ||, where the 
kaon structure is derived in a constituent quark model; 
the model of Ref. || , where some quark and gluon sea is 
considered in addition to the constituent quarks at the 
scale where evolution starts; and most recently, the model 
presented in Ref. , where an expansion of the pion and 
kaon wave- functions in terms of hadron- like Fock states is 
done at the input Q 2 scale. One interesting feature of the 
last one is that the model predicts the parton structure 
of the pion and kaon, at the low Q 2 scale where evolution 
starts, up to a few parameters which must be fixed from 
experimental data. This is a remarkable property since 
leaves the model free of the arbitrariness of the other 
two, where by no means both the form and initial parton 
distribution functions (PDF) are justified. 

In what follows, after a short revision of the model of 
Ref. 0, we shall use it to extract information on the 
pion and kaon parton content from simultaneous fits to 
the available experimental data. 

Following Ref. 0, the n~ and K~ wave functions can 
be written as 













= aj 7T ) 




ba£| 



K 



\K-)+af 



q \k°k~) + 

K-g) + 



K 

■ a 2 



(1) 
■(2) 



at some low Q 2 scale, where we have neglected higher 
order contributions involving heavier mesons and fluctu- 
ations to Fock states containing more than two mesons. 
These fluctuations should be far off-shell and they can be 
safely ignored at this point. The first terms in Eqs. ([!])- 
(^) are the bare meson states, which are formed only by 
valons ||. The following terms are fluctuations whose 
origin can be traced back to the processes v — > v + g and 
g — > q + q followed by the recombination of the perturba- 
tive q q pair with the valons, v, to form the hadron-like 
structure. As discussed in Ref. [0, we assume that in- 
teractions between a quark or antiquark and a valon of 
the same flavor do not form a neutral, unflavored, virtual 
meson structure but annihilates nonperturbatively to a 
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gluon. Additionally, in-meson hadrons are assumed to be 
formed only by valons. 

Thus, these fluctuations are an effective representation 
of the nonperturbative -intrinsic- sea of quarks and glu- 
ons which should provide the necessary binding among 
constituent quarks to form hadrons || |{|. DGLAP evo- 
lution to higher Q 2 generates the perturbative - extrinsic- 
sea of qq pairs and gluons. 

It is worth noting that, on very general basis, individ- 
ual hadrons in the \MM') Fock states are colored. The 
same is also true for the \M g) fluctuations in the sec- 
ond term of Eqs. (|l|)-(§) as far as the gluon is in a color 
octet state. However, the fluctuation itself is colorless. 
Notice that, if the two components of a generic fluctua- 
tion are in the 8 representation of color SU(3), and since 
8(g>8 = 1®8®8® then there is a singlet (colorless) rep- 
resentation where the fluctuation can be accommodated. 

Coefficients af' 1 ; M = ir,K; i = 1,2,3,...; in 
Eqs. (Ey)-(pp are constrained by probability conservation: 

Ei\<\ 2 = i. 

Since valon distributions in pions and kaons are related 

by 
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due to isospin invariance, then at the Q 2 scale, parton 
distribution functions can be written as 

Un(x) = d^ix) = \a%\ 2 v 7r (x) + \al\ 2 P 7rg (g>v 7V 

+ K\ 2 P K K®v k (6) 
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for pions, and 

u K (x) = \a£\ 2 v k (x) + | of \ 2 P Kg <8> v k 

+ |of \ 2 P vK ® «* (9) 

s K {x) = \a£\ 2 v s/k (x) + |af \ 2 P Kg ®v s/k 

+ \a%\ 2 P K Tr®v s/k (10) 

d K {x) = |fl^| 2 P w jc<8» w (11) 
d K (x) = \4\ 2 P Kv ®v k (12) 
g K (x) = |af \ 2 P gK (x) (13) 

for kaons, where 
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The meson probability density Pmm'(x) in the \MM') 
fluctuation has been calculated in Refs. || [To). It is given 
by 

Pmm>(x)= f —F(y,z)R{x,y,z) (15) 

Jo y Jo z 



with 
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In Eqs. v q and q' are the valon and the quark or 

antiquark distributions which will form the meson M in 
the \MM') fluctuation. The q' distribution is generated 
through the gluon emission from a valon followed by the 
q' q' pair creation. Thus its probability density is jllj] 
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where P qg {z) and P gq (z) are the Altarclli-Parisi splitting 
functions p3, 
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The only scale dependence appearing in Eq. ([iS|j arises 
through a st (Q 2 ). Since the valon scale is tipically of the 
order of Q 2 ~ 0.64 GeV 2 then the q 1 q' pair creation 
can be safely evaluated perturbatively because a 2 t /(27r) 2 
is still sufficiently small. The normalization constants 
a, and N in Eqs. @-(^) and @ contribute to 
the global normalization coefficient of the corresponding 
Fock state fluctuation in the expansions of Eqs. (fj])-(@)- 

Momentum conservation also requires 
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a condition which relates the in-meson M and M ' proba- 
bility densities. Additionally, hadrons in the \MM') fluc- 
tuation must be correlated in velocity in order to form a 
bound state. This imply that 
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fixing the exponent a in Eq. ([16]). Notice also that P 9 m 
is calculated from the "recombination" of an antiquark 
with a valon of the same flavor fl . Then, formally, the 
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/ PMM'{y)v q /M ( — ) (14) p M corresponds to the would be tt° distribution in a 
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is the probability density of the nonperturbative contri- 
bution to theparton distribution coming from the \MM') 
fluctuation H liotl. 



hypothetical |M7r ) fluctuation. 

With the 7T and K PDF given in Eqs. (|)-(|) and 
Eqs. (|s|)-([l3|) we proceed to fit simultaneously the ex- 
perimental data on the and uk/u^ by the E615 j^j 
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TABLE I: Coefficients of the pion and kaon PDF obtained 
from simultaneous fits to experimental data from the E615 |^] 
and NA3 pi Experiments. The x 2 /d-o.f. of the fit is 0.882. 

— 2 I K 1 2 

Note that coefficients | | and a 2 are not independent but 
fixed by probability conservation. 
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and NA3 || Collaborations respectively. To that end we 
parametrize [p| 
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where the last two are related one another by momentum 
conservation; and construct a x 2 function as 
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to be minimized. f3{a,b) in Eqs. (|23|)-(|25|) is the Beta 
function, which gives the correct normalization. 

We have a total of 8 free parameters to be fixed by the 

fit, namely, a^' K , a[ ,K and the four exponents in the 

valon distributions of Eqs. (|23|)-(|25|). The fit procedure 
started from an arbitrary set of parameters, for which the 
low Q 2 7r and K PDF were calculated and evolved to the 
Q 2 scale of the E615 and NA3 experiments respectively. 
Then the \ 2 was calculated. This procedure was repeated 
until a minimun of \ 2 was found. 

Results of the fit are shown in Table | and displayed 
in Fig. [l]. As evidenced by the figure, the quality of the 
fit is quite good. 

Once the parameters in the pion and kaon wave- 
functions are known, we build the full set of PDF. They 
are shown for both, the tt~ and the K~ , in Figs. |^ and 
respectively at Q 2 = 30 GeV 2 . To obtain the pion and 
kaon PDF we have used the central values of parameters 
in Table Note that, since |aj| 2 ~ 0, then the quark sea 
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FIG. 1: Fit to experimental data from the E615 H and 
NA3 [[| experiments. Upper: the uuju^ ratio as a function 
of x as measured by NA3. Lower: the distribution as a 
function of x as measured by the E615 Collaboration. 



of the pion is nearly SU(3)-flavor symmetric |L|. How- 
ever, in the quark sea of the kaon the SU(3) symmetry is 
largely broken due to the big contribution of fluctuations 
to the intrinsic sea. Another interesting feature of the 
kaon structure is the d — d asymmetry of the K~ sea. 
This is due to the fact that d is a valon in the 7r while 
d is a valon in the K in the \Ktt) fluctuation. Isospin 
symmetry tell us that (d — d) K (x) = (d — d) K+ (x) — 

(u — u) K (x) = (u — u) K °(x). The u and s sea are, 
however, equal in the K~ . Once again, via the isospin 
symmetry, we get u = s in the K + sea, d — d in the K 
sea and d = s in the K° sea jl3| . 

The PDF of tt + , K+, K° and K° can be obtained just 
by using isospin symmetry. The tt° structure, however, 
deserves separate considerations. In fact, the ir° wave- 
function at the low Q 2 scale can be writen as S 

= 6o|7r )+6 1 |vT%)+6 2 |7r- 7 r + ) 

\K°K°)]+. 



(27) 



Thus, contrarily to charged pions, where the intrinsic sea 
is formed only by gluons and eventually strange quarks, 
in the tt° the intrinsic sea is formed also by uu and dd 
quarks due to the |7r + 7r~) fluctuation. It should be noted 
also that, as the 1 7r° g) and the |7r + 7r~) fluctuations have 
the same origin, namely the splitting of a gluon to a uu or 
dd pair, then it can be assumed that |6i| 2 ~ I&2I 2 , thus 
possibly reducing the intrinsic gluon sea due to proba- 
bility conservation. This indicates a remarkable differ- 
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FIG. 2: Parton distribution functions in charged pions at 30 
GeV 2 . Upper: valence quark distributions (xd(x) = xu(x) = 
xiiv(x) +xu s {x)). Lower: sea quark u n — — s w = s n (full 
line), and gluons (dot line). 



ence between the structure of charged and neutral pions 
within the model. 

We want to stress that we have determine the com- 
plete structure of charged pions and kaons from a min- 
imal set of measurements of the n~ and K~ PDF just 
by extracting the three valon distributions v n , Vk and 



v s / k and the parameters 



ir,K 



and 



from exper- 



imental data. In this sense the model has an interesting 
predictive power. Note that, as a matter of fact, the ex- 
perimental information one can get on the kaon structure 
is only on the light valence quark distribution. The mea- 
surement of the strange and even sea quark distributions 
in kaons is not possible due to practical reasons. Actu- 
ally, strange and sea quarks only contribute to the total 
Drell-Yan dilepton cross section through valence-sea and 
sea-sea qq anihilation. Then their contributions are small 
and cannot be easily separated. A similar difficulty arises 



to measure the sea parton distributions in pions, where, 
once again, valence-sea and sea-sea contributions to the 
Drell-Yan cross section are small than the valence- valence 
ones. 

Moreover, the model predicts the structure of pions 
and kaons at the low Q 2 , scale, where perturbative QCD 
evolution starts. This gives a plausible solution to the 
long standing problem of the origin of the valence-like 
sea quark and gluon distributions nedeed, at the low Q 2 
input scale for evolution, to describe experimental data 
on hadron structure. 
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FIG. 3: Parton distribution functions in kaons at 30 GeV 2 . 
Upper: valence quark distributions Uk = u v + u s (full line) 
and + s s (dashed line). Lower: sea quark and gluon 

distributions Uk = Sk (full line), dk (dashed line), dk (dot 
dashed line) and gluons (dot line). 
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